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ABSTRACT 

Infrared-Faint Radio Sources (IFRS) are objects which are strong at radio wave- 
lengths but undetected in sensitive Spitzer observations at infrared wavelengths. Their 
nature is uncertain and most have not yet been associated with any known astrophys- 
ical object. One possibility is that they are radio pulsars. To test this hypothesis we 
undertook observations of 16 of these sources with the Parkes Radio Telescope. Our 
results limit the radio emission to a pulsed flux density of less than 0.21 mJy (assum- 
ing a 50% duty cycle). This is well below the flux density of the IFRS. We therefore 
conclude that these IFRS are not radio pulsars. 
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1 INTRODUCTION 

Infrared-Faint Radio Sources (IFRS) are objects which typ- 
ically have flux densities of several mJy at 1.4 GHz, but are 
undetected at 3.6 /im using sensiti ve Spitzer Space Te lescope 
observations with /ijy sensitivities l|Norris et al ]|2006l '). They 
were an unexpected discovery in the Australia Telescope 
Large Area Survey (ATLAS), which, covering seven square 
degrees at 1400 MHz to ~30 /iJy, is the widest deep field ra- 
dio survey attempted thus far (Norris et al. 2006; Middelberg 
et al. 2008a). There are ~ 50 IFRS in ATLAS, accounting 
for 2.5% of the current ATLAS source catalogue. Most of 
the IFRS are unresolved at 1.4 GHz and have fiux densities 
of up to tens of mJy. 

The ATLAS project surveys two regions of sky. 
The first region includes the C handra Deep Field South 
(GDFS) |g laccom et all 120011 1 which encompasses the 
Great Observatories O rigins Deep Survey (GOODS) field 
l|Giavalisco et al.] bOOll ) . The second region coincides with 
the European Large Area ISO Survey-Southl (ELAIS- 
Sl). Both regions have also been covered by the Spitzer 
Wide- Area Infrared E xtragalactic (SWIRE) survey program 
iLonsdale et all 120031 '). The second region is close to (but 
does not overlap) the region covered by the recent Parkes 
High Latitude Pulsar Survey ( Burgav et al.,. 20061 ). 

The primary goals of ATLAS are to trace the cosmic 
evolution of active galactic nuclei (AGN) and star-forming 
galaxies. The detected sources are expected to follow rela- 
tively well-known spectral energy distributions (SED) that 



predict infrared emission detectable by Spitzer from most 
strong radio sources. Since the discovery of IFRS in 2006, 
several publications have attempted to understand their na- 
ture and emission mechanisms. Very long baseline inter- 
ferometry (V LBI) observations of a total of six IFRS by 
iNorris et al] ((2003) and Middelberg et al. (2008b) have re- 
sulted in the detection of high-brightness t emperature cores 
in two IFRS. Recentlv. iHuvnh et al.l l|2010l ) investigated the 
four IFRS in the GOODS field, for which ultra-deep Spitzer 
imaging had recently become available, and attempted to fit 
template SEDs corresponding to known classes of galaxies 
and quasars. Norris et al. (2010) have extended this work 
using deep Spitzer imaging data from the Spitzer Extra- 
galactic Representative Volume Survey project (Lacy et al. 
2010), and have shown that most IFRS sources have ex- 
treme values of S20/S3.6, which can be fitted by a high- 
redshift radio-loud galaxy or quasar. They also stacked the 
deep Spitzer data to show that a typical IFRS must have a 
median 3.6 /^m fiux density of no more than 0.2 /iJy, giv- 
ing extreme values of the radio-infrared fiux density ratio. 
Middelberg et al. (2011) have found the spectra of IFRS are 
remarkably steep, with a median spectral index of —1.4 and 
a prominent lack of spectral indices fiatter than —0.7. We 
note that the systematic steepening of spectral indices can 
not be the result of variability. 

While unidentified radio sources have been well-known 
for many decades, they are usually identified once suffi- 
ciently sensitive optical/IR observations become available. 
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The IFRS differ from previous classes of unidentified sources 
in their extreme ratio of radio/IR flux density. For example, 
the ratio of 20 cm to 3.6 /im flux density of the IFRS stud- 
ied here is typically 1000, and is 8000 in the case of ES247 
(Norris et al. 2010). This may be compared with the ratio 
of 10 expected for starburst galaxies and 100 for a typical 
radio-loud quasar template l|Elvis et al.l 1 19941 '). They bear 
so me similarities to the hiRh -redshift radi o galaxies studied 
bv lSevmour et aL I l|2007h andQ arvis et al.l 12009 ) but appear 
in many cases to be even more extreme. Hence, the nature of 
IFRS remains unclear, largely because they are invisible to 
even the most sensitive optical/infrared observations. While 
at least some are likely to be high-redshift radio galaxies 
or quasars, it is also likely that the class may encompass 
several types of object. Other candidates for IFRS include 
heavily obscur ed AGN, str ay radio lobes, radio relics and 
radio pulsars (|Norris et al.l 12006. 2007 ; Gar n & Alexander 
l2008l : Middelberg et al. 2008b'). iNorris et all |2007|) propose 
that IFRS may even represent a phase of AGN evolution. 

For 16 of these sources the measured flux densities and 
spectral indices are consistent w ith those measur ed for pul- 
sars in our own galaxy (Manchester et al] 120051 1'] . As the 
ATLAS survey regions are well away from the Galactic 
plane, it is unlikely that all 16 of these sources are pulsars, 
but the possibility remains that some are. If so, then it is 
important to identify them so that (a) they can be removed 
from the extragalactic source statistics, and (b) so they can 
be accounted for in models of pulsar populations. 

However, it is difficult to estimate how many pulsars 
would be expected in any given region of the sky. The 
most recent pulsar s urveys have covered reg ions close to the 
Galactic plane (e.g.. [Manchester et al.ir2001i ). These pulsar 
surveys have discovered a large number of new pulsars, but 
relatively few millisecond pulsars. The Galaxy is likely to 
contain a similar number of normal and millisecond pulsars 
l|Lvne et al.lll99j ). but the millisecond pulsars are likely to 
have travelled significant distances from the Galactic plane 
and are difficult to find due to their relatively low lumi- 
nosities, their fast spin-rates and because they are com- 
monly in binary systems, which disrupts the regular peri- 
odicity of their emitted pulses. The only deep 20-cm pul- 
sar survey covering a large area that is far from the Galac- 
tic plane is the Parkes High Latitude Pulsar Survey, which 
covers Gal actic longitudes 22 0° < I < 260° and latitudes 
|6| < 60° l|Burgav et al.ll2006l ). with 6456 pointings of 265 
seconds each. This survey detected 42 pulsars with a survey 
sensitivity of pulsed flux density ~ 0.5 mjy. Even though 
this survey made relatively few discoveries, the pulsars de- 
tected away from the Galactic plane were of great astrophys- 
ical interest. They included the first double pulsar system 
PSR J0737-303 9 (iLvne et al.ll2004) and three other mil- 
lisecond pulsars. iBurgav et al.l (|2006l i discussed the Galactic 
latitude distribution of their detections in detail, and sug- 
gested a roughly uniform distribution. However, this survey 
was relatively insensitive to millisecond pulsars, because of 
the restricted sampling time and number of frequency chan- 
nels, and so the total number of millisecond pulsars detected 
was small. Hence, it is currently impossible to make a rea- 
sonable estimate of the number of pulsars likely to be de- 
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tected per square degree at high Galactic latitude. With a 
survey sensitive to millisecond pulsars it seems reasonable 
that at least one detectable pulsar exists in the ATLAS sur- 
vey region of roughly seven square degrees, but because this 
number is so poorly known, it adds a further motivation for 
this study. 

It should be noted that detections and discoveries of 
pulsars in radio continuum data are not without precedent. 
Exa mples include PSR B19374-21, the first millisecond pul- 
sar (|Backer et al.lll982l ). PSR B1821-24, found during a 
VLA imaging search of globular clusters (|Hamilton et al.l 
Il985l ). PSR B1951-I-32, initially identified as a steep- 
spectrum polarise d point sou rce central to the supernova 
remnant CTB 80 (|Stromlll987i l. and PSR J0218+4232, orig- 
inally uncovere d by the Westerbork Synthesis Radio Tele- 
scope (WSRT) jNavarro et al.ll 19951 ). 



2 OBSERVATIONS 

We selected 16 of the IFRS based on their flux density 
and spectral index. The sources range in flux density from 
0.5 mJy to 23 mJy and possess spectral indices from —2.4 
to —0.2. An explanation of the calculation and derivation 
of these values, including the problems associated with res- 
olution and scintillation, may be found in Middelberg et al. 
(2011), from which the IFRS values were taken. These pa- 
rameters are certainly similar to those of the known pulsar 
population, which have flux densities that range from 0.01 to 
1100 mJy at an observing frequency close to 1400 MHz and 
have known spectral indices ranging between —3.5 and +0.9 
with a median of —1.7 (Manchester et al..i,2005y It should be 
noted that the lower value of the pulsar flux density range 
is not due to any physical characteristic of pulsars, but is 
instead due to the physical limitations in the sensitivity of 
our instrumentation. 

In Table[T]we tabulate the source names, positions, flux 
density and spectral index of the IFRS. For each source we 
carried out a 35 minute observation using the Parkes 64-m 
radio telescope equipped with the 20-cm 13-multibeam re- 
ceiver (although only the central beam was utilised for our 
analysis), with 340 MHz of bandwidth centred on 1352 MHz. 
We used the Berkeley Parkes Swinburne Recorder (BPSR), 
a high resolution digital filterbank, providing 870 frequency 
channels, 2-bit samples every 64 /is. This observing setup 
is almost identical to t he High Time Re solution Universe 
Pulsar Survey project (|Keith et al.l |2010| ). using a similar 
bandwidth, observing frequency backend system and sam- 
pling rate, with high sensitivity to millisecond pulsars. The 
data were recorded to magnetic tape for off-line processing. 
The data were p rocessed using the Hitrun pipeline (see 
iKeith et al.ll201ol for details) to search for periodic signals in 
dedispersed time series with trial dispersion measures (DMs) 
in the range of to 1000 cm~^ pc. 

We determine the theoretical sensitivity of our observa- 
tions using the radiometer equation, which gives the funda- 
mental limiting fiux density of the central beam. 
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In this equation, a is the cutoff S/N for a positive detec- 
tion, Tsys is the system noise temperature, Tsky is the sky 
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Table 1. Parameters of the infrared- faint radio sources selected for analysis as potential pulsars. All flux densit y and sp ectral index data 
is obtained from Middelberg et al. (2011), with the exception of ES1259 (Middelberg et al. 2008a) and CS255 l|Norris ct al. 2006. 1. 
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noise temperature, G is the system gain, B is the observing 
bandwidth. Tabs is the integration time of the observations 
and W is the fractional pulse width. The typical value of 
W for pulsars is approximately 10 percent of the pulse pe- 
riod, although it can be as high as 5 perc ent. Scaling our 
parameters from those of iKeith et al.l (j2010l ). setting cr = 8, 
Tsys = 23 K, Tsky = 0.7 K and G = 0.735, we arrive at 
a sensitivity to pulsed emission of Smin = 0.07 mJy for a 
10 percent pulse duty cycle, or Smin — 0.21 mJy for a 50 
percent pulse duty cycle. As all of our chosen IFRS candi- 
dates possess continuum flux densities above the higher of 
the two flux density limits, if any of them are indeed pulsars 
they should be easily detectable by our observations at the 
8cr confidence level. 



3 RESULTS AND DISCUSSION 

The folded pulse profiles (intensity as a function of pulse 
phase) were viewed by eye for all candidates with a signal- 
to-noise ratio above 8a. The observed pulse profiles for each 
candidate and their pulse periods were clearly caused by 
radio frequency interference or processing artefacts. For ex- 
ample, one candidate with a signal-to-noise ratio of 13a was 
detected with a frequency of 50.01 Hz and a sinusoidal pulse 
profile, indicating RFI originating with the mains power sup- 
ply. We note that the minimum fiux density of the chosen 
IFRS is 0.5 mJy, and the median fiux density of the IFRS 
is approximately 4.4 mJy. If these fiux densities are due to 
pulsars, then with a duty cycle of W = 0.5 and an 8a de- 
tection level of 0.21 mJy, the faintest of the IFRS should 
have been detected at a confidence level of 19a with most of 
the sources ruled out by more than 100a. We conclude that 
none of the IFRS observed in this project show emission 
consistent with being a radio pulsar. 

The absence of pulsar detections in this project empha- 
sises the difficulty of identifying pulsar candidates through 
continuum surveys. Despite the four detections described in 
Section[TJ the vast majority of continuum (i.e., low time res- 



olution) pulsar surv eys have found no positive results (e.g. 
iKaplan et al]|2000l and Crawford et al. 2000). Also, while 
our analysis attempted to take advantage of the typically 
steep values of pulsar spectral indices in selecting our tar- 
gets, the principle means of searching for pulsar candidates 
in continuum surveys has been to identify sources with sig- 
nificant linear a nd/or circular polarisation (Han ct al. 1998|; 
iHan et al.|[2009h . Further studies using these selection crite- 
ria may yield better results. 

Finally, the Evolutionary Map of the Universe (EMU) 
survejQ, to be carried out using the Australian Square Kilo- 
metre Array Pathfinder (ASKAP), is likely to detect 70 mil- 
lion radio sources with a sensitivity limit of ~ 10/iJy. Over a 
million of these sources are likely to be currently undetected 
IFRS. However, at this point in time, we cannot make a pre- 
diction as to the number of these IFRS which may be pul- 
sars, or indeed the fraction of pulsars expected in the total 
EMU source count. These values, as well as the techniques 
which may be used to more accurately isolate candidate pul- 
sars from large-scale continuum surveys such as EMU, are 
the subject of ongoing research and are to be published in 
later work. 



4 CONCLUSIONS 

We have searched for short-term radio pulsations originating 
from 16 of the ATLAS infrared faint sources and find that 
pulsed emission cannot account for the observed radio fiux 
density. Since the sensitivity limit of our observations is well 
below the observed flux densities of the chosen IFRS, it is 
unlikely that any of these enigmatic sources are simply close- 
by pulsars in our own Galaxy. Hence, the nature of these 
sources is yet to be determined. 



^ http:/ /www. atnf.csiro.au/people/rnorris/emu/ 
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